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ABSTRACT 

Cygnus 0B2 is the nearest example of a massive star forming region, contain- 
ing over 50 0-type stars and hundreds of B-type stars. We have analysed the 
properties of young stars in two fields in Cyg 0B2 using the recently published 
deep catalogue of Chandra X-ray point sources with complementary optical and 
near-IR photometry. Our sample is complete to ~1 (excluding A and B-type 
stars that do not emit X-rays), making this the deepest study of the stella r prop - 



erties and star formation history in Cyg 0B2 to date. From ISiess et al.l (120001 ) 
isochrone fits to the near-IR color- magnitude diagram, we derive ages of 3.5tiJ)^ 
and 5.251}'q Myrs for sources in the two fields, both with considerable spreads 
around the pre-MS isochrones. The presence of a stellar population somewhat 
older than the present-day 0-type stars, also fits in with the low fraction of 
sources with inner circumstellar disks (as traced by the i^'-band excess) that we 
find to be very low, but appropriate for a population of age ~5 Myrs. We also 
find that the region lacks a population of highly embedded sources that is often 
observed in young star forming regions, suggesting star formation in the vicinity 
has declined. We measure the stellar mass functions in this limit and find a 
power-law slope of F = —1.09 ± 0.13, in good agreement with the global mean 
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value estimated by iKroupal (120021 ) . A steepening of the slope at higher masses 
is observed and suggested as due to the presence of the previous generation of 
stars that have lost their most massive members. Finally, combining our mass 
function and an estimate of the radial density profile of the association suggests a 
total mass of Cyg 0B2 of ~ 3 x 10^ M©, similar to that of many of our Galaxy's 
most massive star forming regions. 

Subject headings: stars: early-type - stars: formation - stars: pre-main sequence 
- Galaxy: open clusters and associations: individual (Cygnus 0B2) 



Introduction 



Stars do not form individually but in groups, with sizes ranging from smal l star forming 
regio ns (SFRs) to the massive superstar clusters seen in merging galaxies ( iLada fc Lada 
20031 ). Our current understanding of star formation is mainly derived from observations of 
nearby small star forming regions such as Taurus or Ophiuchus. However the conditions of 
star formation vary considerably between these regions and the massive SFRs that contain 
hundreds to thousands of OB stars and millions of low-mass stars. High stellar densities, 
strong stellar winds from OB stars and a large UV flux are likely to influence the products 
of the star formation process such as the initial mass function (IMF), the binary fraction 
and the evolution of protostellar disks. Our understanding of the influence of environment 
on these products is currently poor, and needs to improve if we are to develop a global 
understanding of star formation that can relate the small-scale properties with the evolution 
of galaxies and the survival of planetary systems. 

The majority of massive SFRs are found in external galaxies, or close to the center of our 
own Galaxy. Their rarity and great distances therefore make them difficult ob jects to study . 
The one exception to this is Cygnus 0B2, which, at a distance of only 1.45 kpc ( lHansonll2003l ) 
is the clos est known massive S FR. Early surveys of the region revealed the presence of ~300 
OB stars ( jReddish et al.lll967l ). a number that has grown as deeper surveys have penetrated 



the high exti nction in its direction (e.g. iMassey &: Thompson 



1991; Comeron et al. 2002 



HansonI |2003[ ) . There are currently 65 known 0-type stars in the region making it the 
largest known concentration of spectroscopically confirmed 0-type stars. 



From a near-IR photometric study iKnodlsederl (120001 ) estimated the region contained 
over 100 0-type stars, ~2600 OB stars, and had a total mass of 4-10 xlO^ M©. However, 
the complications of background subtraction in the Galactic plane, combined with strongly 
variable extinction in the vicinity of Cyg 0B2, give cause to doubt such an estimate. Further- 
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more the question of the extent of Cyg 0B2 is made more comple x by the growing ev idence 
for stars older than the 1-3 Myr age appropriate to the OB stars: iDrew et al.l ( 120081 ) found 
a concentration of 5- 7 Myr old A stars in 1 sq. deg just to the south of Cyg 0B2, while 
Comeron et al.l (|2008[ ) identified an unclustered spread of evolved stars up to ~10 Myrs old 
over a wider area, 2 degrees in radius. In this study we obtain further evidence that would 
favour the present-day OB stars as the products of the latest but most certainly not the only 
phase of star formation. 

Recent deep photometric surveys of the Galactic plane in the optical and near-IR have 
penetrated much of the extinc tion toward the region and are adva ncing our understanding of 



this massive association (e.g. iDrew et al.l 120081 : IVink et al.ll2008l ). The combination of data 



from these surveys with that at other wavelengths is vital for studying the star formation 
process, particularly in the dense and highly obscured regions of the Galactic plane, and one 
that this work makes particular use of. 

One of the objectives of this work is to measure the IMF down to low masses. The IMF 
is one of the most important measurable quantities in star formation studies. Combined with 
the time-variable star formation rate, the IMF determines the evolution of stellar systems 
ranging from small SFRs to entire galaxies. Mea surements of the IMF have shown it to be 
highly uniform at high and intermediate masses ( Luhman et al.l 2000 : Masseylll998 ) with a 
power law slope above 1 Mq with an index of F = —1.35 ( ISalpeterl 119551) for in termediate 
mass stars, possibly steepening to F = —1.7 for the most massive stars (IScalolll986[ ). At lower 



masses the IMF is less well constrained, but appears to flatten below 0.5 M(7^ and e xhibits a 
'turnover' around 0.1 - 0.3 Mq, with less stars of the lowest masses (lKroupall200ll ). 



Deviations from this uniformity are becoming increasingly rare as deeper and more com- 
plete observations overturn previous measurements and the majority of remaining observa- 
tions have been shown t o be within the statistical uncertainties of the seemingly universal 
mean (e.g. lKroupall2002[ ). While some theoretical models predict a shallower IMF in regions 



of high dens ity or pressure, on ly a small number of observa tions in the vicinity of the Galactic 



Center (e.g. IStolte et al.ll2006[ ) or in starburst regions (e.g. ISmith fc Gallagherll200ll ) support 



this and more measurements of the IMF in massive SFRs ar e necessary to r esolve this. The 
most recent measurement of the IMF in Cyg 0B2, due to iKiminki et al.l ( 120071 ). suggests 
that at high masses it is in fact steeper than the canonical value (though they estimate their 
spectroscopic sample becomes incomplete at masses below ~15 Mq). This further supports 
the need for a thorough and complete measurement of the IMF in Cyg 0B2 down to lower 



masses. 



In this paper we present an analysis and discussion of the stell ar properties of sources 
in Cygnus 0B2 taken from the X-ray derived catalogue presented by IWright fc Drakd ( 120091 . 
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hereafter Paper 1). In Section 2 we discuss the source catalogue and the various cuts apphed 
to it to remove foreground and background contaminants. In Section 3 we present the near- 
IR properties of these sources, in Section 4 we study their X-ray properties and in Section 5 
we derive mass functions for the sources. Finally, in Section 6 we discuss the implications of 
our findings on our understanding of Cygnus 0B2 and compare it to other massive SFRs. 



2. Cyg OB2 membership and foreground contamination 



The location of Cygnus 0B2 in the Galactic Plane introduces the non-trivial issue of 
contamination from foreground and background sources. This problem is further compli- 
cated by the high and variable extinction toward and throughout the Cygnus region. There 
are many methods that have been used to select members of a star forming region while 
minimizing the contamination from the field population, but one of the most effective and 
efficient is to use X-ray emission as a tracer of youth, since you ng pre-MS stars are orders o f 
magnitude more luminous in X-rays than their MS equivalents ( iPreibisch fc Feigelsorul2005l ). 
In Paper 1 a catalogue of 1696 sources X-ray sources was presented, extracted from Chandra 
observations of two fields in the center of Cyg 0B2. The catalogue also in cluded optical an d 
near-IR associations taken from the I PHAS (INT Photomet ric Ha Survey. iDrew et al.ll2005l ). 
2MASS (Two Micron All Skv Survev. Iskrutskie et aD l2006h and UKIDSS (United Kingdom 
Infrared Deep Sky Survey, iLucas et al.ll2008l ) photometric catalogues. This work is based on 
this catalogue and principally the 1501 sources with either optical or near-IR counterparts 
that represents the main stellar sample. 

The completeness of this catalogue is dependent on both the X-ray detections and the 
associations at other wavelengths. In Section H^l we show, by comparison with deeper X-ray 
stellar catalogues, that the completeness limit for this catalogue is ~1 Mq. In Paper 1 it 
was shown that the near-IR photometry combined from 2MASS and UKIDSS observations 
will allow pre-MS stars down to 0.2 Mq to be detected at an extinction of Ay ~ 7. In 
the north-western field, where some UKIDSS data is lacking, it was estimated in Paper 1 
that the 2MASS observations reach a pre-MS depth of 1 M© at Ay ~ 7, though the lower 
extinction in this field means the observations will be slightly deeper than this. 

Since not all young stars emit X-rays we must also be aware of the limits of X-ray 
emission mechanisms. O and early B-type stars (earlier than type B2, >10 M© are expecte d 



Berghoefer et al. 



19971) 



to generate X-rays in their radiatively-driven stellar winds (e.g. 
while, at ages of ~5 Myr, lower mass stars (<3.3 Mq) have substant ial convective envelope s 
that produce X-ray emission through magnetic dynamo activity (e.g. lPallavicini et al.lll98ll ). 



Intermediate mass stars do not have deep convective layers or strong stellar winds and are 
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Fig. 1. — IPHAS (r'— Ha, r' — i') color-color diagrams for all IPHAS sources in the Chandra 
observational footprints (left-hand panels) and only for sources associated with Chandra 
sources (right-hand panels) for the central (upper panels) and north-western fields (lower 
panels). Only sources with r' < 20 and magnitude er rors < 0.1 are shown. Blue dashed lines 
show unreddened main-sequences (jPrew et al.ll2005l ). while red dashed lines are the main- 
sequence tracks for reddenings of Ay = 4.0 and 7.0 (central field) and Ay = 2.5 and 5.0 
(north-western field), illustrating the region approximately encompassing the Cygnus 0B2 
population. Sources removed from the catalogue as foreground sources are shown as red 
crosses. 
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therefore not expected to emit X-rays. Studies that have observed X-ray emission from 
late B and A-typ e stars have attribu ted this to coronal emission from unresolved late-type 
companions (e.g. IStelzer et al.ll2006[ ). The intermediate-mass completeness regime can be 
determined by counting the number of spectroscopically known OB stars in our field of view 
with X-ray detections. We find that the fraction of OB stars detected in X-rays decreases 
from 100% for all 0-type stars to 64% at BO, 36% at Bl and < 20% for all later intermediate- 
mass stars. This suggests we are incomplete for the majority of B-type stars, and also most 
early A-type stars. 

Finally, crowding is likely to have only a minor effect on completeness, because only 
the very central regions of the association are particularly crowded and these were observed 
close to on-axis with Chandra and have complete UKIDSS coverage, therefore allowing them 
the maximum possible spatial resolution in both X-rays and in the near-IR. 



2.1. Removal of foreground sources 



While X-ray observations are excellent at separating young sources from older field 
stars, this method is not infallible since late- type main-sequence (MS) stars of all ages have 



been observed to have X-ray emission (e.g. IVaiana et al. 



198ll ). albeit at generally lower 



levels. Nearby MS stars, although intrinsically fainter, could therefore have similar X-ray 
brightnesses as distant pre-MS stars, and thus be a potential contaminant in our catalogue. 
To remove any foreground sources, we employ one of the interesting features of the IPHAS 
{r'—Ha, r' — i') color-color diagram (CCD) where the main sequence does not redden onto 
i tself like it does in the near-IR CCD, but sweeps out an area in the color-color plane 
(IDrew et al.l l2005l ) . This allows differently reddened stellar populations along the hne of 



sight to be separated (e.g. iDrew et al.ll2008l : ISale et al.ll2009l ). 



In sightlines towards the Cygnus region, iMcCuskey &: HoukI (jl97ll ) found a high con- 
centration of stellar sources between 200 and 300 pc and a strong decline in stellar density 
between 500 pc and 1 kpc, with the stellar populations making up Cygnus 0B2 and Cygnus X 
beyond that. Figure [1] shows IPHAS CCDs for all IPHAS sources and for just those with 
Chandra X-ray counterparts in the are as covered by the Chandra observational footprints. 
It is clear from the CCD, as noted by iDrew et al.l (120081 ). that the majority of foreground 
stars in the direction of Cyg 0B2 are nearby and only mildly reddened. The main Cyg 0B2 
population, as traced by Chandra, can clearly be seen at between Ay = 4-7 in the central 
field and Ay = 2.5-5 in the north-western field. 



It is possible that pre-MS members of Cyg 0B2 with Ha emission may appear above the 
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reddened main-sequences and closer to the foreground population (Ha emission, frequently 
observed in young stars, causes the r'— Ha color to increase). To preclude the removal of 
genuine members of Cyg 0B2 we compare the positions of candidate foreground stars in 
the near-IR color-magnitude diagram (CMD) with a foreground MS. Those sources whose 
positions in both diagrams were appropriate for pre-MS sources at the distance of Cyg 0B2 
and potentially experiencing Ha emission were retained. It should be noted that the majority 
of removed sources fall on a thin strip in the IPHAS CCD that would exhibit a larger spread 
if their positions in the CCD were due to the effects of Ha emission. 

For sources with only near-IR associations this method was not available to us and we 
had to rely on the positions of sources in the near-IR CMD. Due to the smaller effects of 
extinction in the near-IR and the fact that near-IR isochrones redden onto themselves it is 
not simple to separate the foreground population in a similar manner. Instead we may only 
identify sources with very blue colors as likely foreground sources and we have removed a 
small number of outlying sources in this way. In total we removed 46 sources (~3%) from 
our catalogue as likely X-ray emitting foreground sources, which are listed in Table [1] These 
were predominantly removed from the 757 sources (~50% of the original catalogue) with 
IPHAS photometry. Based on the fraction of IPHAS sources removed and the number of 
sources without IPHAS photometry we estimate that a small fraction <1% of foreground 
contaminants may remain in the catalogue, but will be impossible to identify without spectra 
or further photometry. Finally we note that the spatial distribution of sources removed is 
uniform and does not exhibit a clustered distribution as the overall X-ray source distribution 
does (see Paper 1). We therefore confidently remove these sources from further consideration. 



3. Near-IR stellar properties 

We now consider the near-IR properties of the remaining 1497 Chandra sources with 
2MASS or UKIDSS associations that we believe to be members of Cygnus 0B2. We will 
consider only those 2MASS sources with valid quality detections (photometric quality flags 
"A" to "D") and both 2MASS and UKIDSS sources with magnitude errors <1%, in the 
relevant bands. Applying these restrictions to all three bands we are left with 1067 and 
290 sources in the two fields, losing 50 (4%) and 90 (24%) sources in the central and north- 
western fields respectively. We lose a larger fraction of sources in the north-western field 
because the UKIDSS survey has yet to take data in this region and the fraction of 2MASS 
sources in this area with bad quality Kg-henad photometric flags is high. However this should 
not affect the completeness limit of 1 Mq since the available 2MASS photometry is deep 
enough to detect a one solar mass star at the distance and extinction of Cyg 0B2 (c.f. 
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Paper 1). For ease of comparison between the two near-IR surveys and previous work in 
the hterature all UKIDSS photometry w as converted onto th e 2MASS photometric system 
in Paper 1 using the conversions given by iHewett et al.l (120061 ) . These conversio ns were used 



i nstea d of the more accurate extinction-dependent transformations presented by lLucas et al. 



( I2OO8I ) because of the inhomogeneous extinction across the region (the expected systematic 
difference in the use of these transformations will be ~0.03 mag at J, but only ~0.01 mag 
at H and K). 



3.1. Near-IR color-magnitude diagram 



Figure [2] shows (J, J — H) CMDs for sources in the two fields. We use the J — H 
color instead of the H — K color because it is less affected by circumstellar dust and has 
a smaller degeneracy between color and magnitude for pre-MS isochrones. Applying the 
photometric criteria discussed above to the J- and i7-band photometry, there are 1076 and 
325 valid sources, 96% and 86% of sources in each field. Both CMDs show a clear sequence 
with large numbers of sources down to J ~ 19. The greater depth in the central field is due 
primarily to the longer exposure of the Chandra observation, although the completeness of 
the UKIDSS survey in this field is also a benefit. The sources in the central field are also 
redder, as expected considering the greater extinction in this direction. 

To determine stellar properties of the population of X-ray sources we have attempted 



to fit a number of pre-MS isochrones to our data, as shown in Figure [21 We used ISiess et al. 



(I2OOOI ) isochrones for la te spectral types (M < SM©) co nverted onto the 2MASS system 



using the data given by iKenyon fc Hartmanru ( 



we u sed the MS absolu te ma gnitudes in ICramerl ( 



1995) and Besselll (2005 ). For e arlier types 



from iKoornneefl (119831 ) and iKenyon fc HartmannI (119951 ). Goodness-of-fit tests, weighting 



and iMartins et al.l (120051 ) and colors 



each source using its color and magnitude errors, were performed to estimate the best-fitting 
isochrone and mean extinction of the sources in each field. Fits were attempted at a range 
of ages (at intervals of 0.25 Myrs) and reddenings (at intervals in Ay of 0.25). To assess 
the validity of these results, 95% confidence intervals were calculated for all fitted ages and 
extinctions using a bootstrapping technique in which our dataset was resampled 1000 times 
and refit. The range of fitted values provides a reliable estimate of the confidence intervals 
on each fit. Both the age and extinction were left as free parameters during the refitting. 

We note first that, regardless of age, the lower mass stars are better fit by a slightly 
redder isochrone than the higher mass stars. This effect is most prominent for stars < 2 M0, 
and is apparent in both fields, although it is stronger in the north-western field. This 
could be attributed to inaccuracies in the isochrones or near-IR colors for massive stars 
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Fig. 2. — Near-IR CMDs for Cygnus 0B2 sources in the central (left) and north-western 
(right) fields for all s ources that meet the photometric requirements described above. Pre- 
MS isochrones (from lSiess et al.l 120001 ) are shown using the best fit parameters described in 
the text (3.5 Myrs, Ay = 7.5 for the central field and 5. 25 Myis^ Ay = 5.5 for the north- 
western field), and at a distance of 1.4 5 kpc (DM = 10.8, 



Hanson 



2003). For masses above 



8 Mq we us e absolute magnit udes from Martins et al. ( 2005h an d Cramer ( 1997 ) and near-IR 
colors from iKoornneej (1983) and Kenyon Sz Hartmann (1995), for the O- and B-type stars 



respectively. The MS and pre-MS were c o nvert ed on to the 2MASS system using the data and 
relations given by iKenyon fc HartmannI (119951) and iBesselll (120051 ) and then reddened using 
a s tandard Galac t ic ext incti on law (Howarthl Il983l ) and relative extinction relations given 
by ISchlegel et al.l (119981 ) and iHansonl (120031 ). Mean photometric errors for each magnitude 
interval are shown as red crossed along the right-hand side of each figure. 
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(particularly since a set of empirical 2MASS colours for massive stars does not exist in the 
literature), but we note that such a discrepancy has not been noted for other star forming 
regions. The use of two different sources for the near-IR photometry raises the question 
of whether the different photometric systems may have caused this discrepancy, but since 
the discrepancy is apparent in both fields and the north-western field uses primarily 2MASS 
photometry, the different photometric catalogues are unlikely to be re sponsible. The effect 
is also visible in the CMD presented by lAlbacete Colombo et al.l (120071 their Figure 4), who 
used only 2MASS photometry in their analysis of the central Chandra field. The origin of 
this discrepancy is unclear. The high and variable extinction toward Cyg 0B2 and complex 
star formation history (see discussion in Section 16. ip makes this difficult to resolve. Despite 
the apparent difference in reddening however, we show below that the difference between the 
fitted extinctions of the high- and low-mass populations are not significant. 

Because of this we fit the high- (8-20 M0) and low-mass (<8 M©) populations separately. 
The extinction was allowed to vary for both fits, while the age of the region was only varied 
for the low-mass population because the high-mass population should already be on the 
MS and therefore, excluding the highest mass stars which may have evolved off the MS 
and which we do not include in our fit, are insensitive to age variations. In the central 
field the high-mass stars are best fit at a reddening of Ay ~ 7.2b'^^, with the lower-mass 
stars redder with Ay ~ T.SIq^j. In the north-western field the high- mass stars are fit with 



Ay ~ 4.75_q;75 with the low-mass stars having Ay ~ 5.5_]^;o. Since there is no degeneracy 
between extinction and assumed distance, these values are independent of our choice of 
distance modulus. Compared to the mean extinctions calculated in Paper 1 from X-ray 
spectral fits these values are slightly larger, though within the uncertainties calculated here. 

When fitting the low-mass compon ent with pre-MS isochrones we attempted fits at both 
the older distance of 1.8 kpc fDM=11. 2.lMassev fc Thompsorull99ll ) and the revised distance 
of 1.45 kpc (DM=10.8) estimated bv iHansonl J2003h . Using DM=11.2 the low-mass stars 
are best fit with a 2.5j]]^'25 pre-MS isochrone, in agreement with pre vious studies based 



on the evolutionary state of the most massive stars in Cyg 0B2 (e.g. iHansonl |2003| ). At 



DM=10.8 the population is better fit with a slightly older 3.5l?'n^ Myr iso c hrone (as shown 



in Figure [2]). This is older than the upper limit determined by Hanson ( 20031 ) based 
main-sequence 06 and 07 dwarfs, but within the 95% confidence interval. 



on 



In the north-western field we find consistently older ages for the stellar populations. At 
DM=11.2 the low-mass population is best fit at an age of 4.0l}'25 Myrs and at DM =10.8 at 
an ag e of 5.251};q Myrs. For both fields the use of the closer distance proposed by [Hanson 
(I2OO3I ) requires an older pre-MS isochrone to fit the data. These ages will be discussed in 
more detail in Section 16.11 
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3.2. Near-IR color-color diagram 



Figure |3] sliows tlie {J — H, H — Kg) CCD for tlie 1357 sources tliat meet our plio- 
tometric requirements (fewer stars tlian present in the CMD in Figure [2l on account of 
requiring satisfactory K., rn agnitudes also). For reference we also show the M S track from 



Kenvon fc HartmannI (119951 ) and the classical T-Tauri star (CTTS) locus from iMeyer et al. 



( 1l997l ). both with reddening vectors extending from them. As noted above the majority of 
sources trace a reddened main-sequence, with a number of sources significantly more red- 
dened still. The higher photometric errors of the 2MASS-dominated north-western field 
photometry is also evident when comparing the two fields and from the mean color error 
shown in each figure. In the central field there are two sources with large J — H colors and 
low H — Kg colors that do not correspond to plausible stellar colors and as such may have 
spurious photometry. 

The near-IR CCD is often useful for estimating the evolutionary class of young stellar 
objects because sources with circumstellar material will be shifted to higher if — colors 
due to near-IR excesses in the Kg band. These are Class ll objects ( iLadal 119871 ) . such as 
T-Tauri or Herbig Ae/Be stars, which are pre-MS stars with circumstellar accretion disks 
that are less evolved than the Class ill objects that have lost their inner circumstellar disks 
and clearly dominate our sample. We find 63 and 23 sources that fall within this region of 
the color-color plane, 5.9% and 7.9% of our catalogue for the central and north-western fields 
res pectively. The fraction f o r the central field is slightly larger than the 4.4% determined 



by lAlbacete Colombo et al.l (120071 ). which can likely be attributed to our use of UKIDSS 
data to probe lower- niass stars that retain circu mstellar disks for longer periods of time 
(e.g. iHillenbrand et al.lll998t iHaisch et al.ll2001al ). To the right of the reddened CTTS is a 



region typically occupied by less evolved Class I protostars, or evolved stars, with significant 
amounts of circumstellar material. We find only one object in this region, and from its 
isolation and lack of IPHAS photometry we are unable to determine whether it is a true 
member of Cyg 0B2 or a background source. 

The fractions of Class I and ll sources are significantly lower than that found from studies 



of oth er similar age star forming regions in the near-IR such as M17 ( 37%. iHoffmeister et al. 



2008h and the Orion Nebula Cluster (>55%. lHillenbrand et a0l998[ fl. lAlbacete Colombo et al 
( 120071 ) suggested that the low fraction of circumstellar disks in Cyg 0B2 was due to the strong 
UV radiation field from the large number of OB stars in the vicinity that causes photoe- 



^We note that while different studies often use different criteria for defining near-IR excesses and quote 
excess fractions over different stellar mass ranges, the clearly higher fraction of excess sources in these regions 
compared to Cyg 0B2 is worthy of attention and explanation 
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Fig. 3. — Near-IR CCDs for Cygnus 0B2 sources in the central (left) and north-western 
(right) fields for all sources that meet the photometric requirements described above. The 
Kenyon fc HartmannI (119951 ) main-sequence is shown as a blue l ine, t ransformed onto the 
2MASS photometric system using the relations given by iBesselll (120051 ). Reddening vectors 
with length Ay = 15 mag are shown extending from the MS at types K7V (M ~ 0.6 M©) 
and M6V (M ~ 0.2 Mfn) using the redd ening relations described in Figure [21 The classical 
T-Tauri star locus of iMeyer et al.l (119971 ) is shown as a straight red line, also with a reddening 
vector extending from it, and sources falling within this region are marked as crosses. The 
highly reddened object in the central field is shown as a large dot. The mean photometric 
errors on the two colors is shown in the top left corner of each figure in red. 
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vaporation of the disk material. Although predicted by theory (e.g. I Johnstone et al.lll998l ). 
an inverse relationship between disk fractions and UV radiation field has only been rn ildly 



observed in star forming regions to date (e.g. iBalog et al.l 120071 : iGuarcello et al.l 120071 ). To 



a certain extent our results support such a theory, with the fraction of circumstellar disks 
lower in the center of the star forming region where the UV radiation field from OB stars 
is stronger. However the difference is both statistically small and this ir iterpretation over- 
simplifies the dynamical and 3-dimensional structure of the region (e.g. iBalog et al.l |2007| ) 
and therefore does not allow an accurate comparison between disk fractions and the UV radi- 
ation field. However, If Cyg 0B2 contains a significant fraction of sources with ages > 5 Myr , 
then the disk fractions of 5 - 10% found here are perfectly reasonable (e.g. lHillenbrandll2005l ). 
Given the growing evidence for an older generation of star formation in the region (see dis- 
cussion in Section 16. ip we suggest that this is a more likely explanation for the low fraction 
of circumstellar disks than invoking UV photoevaporation from OB stars. 

The use of near-IR colors is however an imprecise method of tracing the population 
of disk bearing stars since a i^^-band excess will only reveal hot dust from the inner parts 
of a circumstellar disk. The combination of near- and mid-IR photometry is necessary to 
produce a more complete picture of the evolution of circumstellar disks in a star forming 
region. Forthcoming observations of the Cygnus region from the Spitzer Space Telescope will 
allow this matter to be probed in more detail. 



4. X-ray stellar properties 

In this section we discuss how the integrated X-ray properties of sources in Cyg 0B2 
may be used to study the star forming region as a whole, and make comparisons with the 
properties of X-ray sources in other star forming regions. 



4.1. Median X-ray photon energy 



Feigelson et al.l ( l2005l ) showed that the median energy of X-ray sources is a reliable indi- 



cator for the absorbing column density, log Nh, obtained from spectral fitting procedures. In 
FigureHlwe show the distribution of median energies of Chandra sources in Cy gnus 0B2 com- 
pared to that of the well-studied Chandra Ultra-deep Orion Project (COUP, iGetman et al. 



2005bl ) of the Orion Nebula Cluster (ONC). The peak median energy of Cygnus 0B2 sources 



is higher than that of Orion because of the higher extinction toward Cyg 0B2 compared to 
Orion. 
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The COUP sample also shows a smaller population with a higher median energy around 
3 keV that is indicative of a population of embedded sources with higher extinction. However, 
such a second population is not clearly visible in the Cyg 0B2 distribution with the high 
median energy sources appearing to constitute the tail-end of the standard distribution. 
This can be quantified as the fraction of sources in the two catalogues that would be classed 
as "heavily obscured" by some authors, i.e. with {E) > 3.0 keV. In Cyg 0B2 only 10% 
of the 1696 sources fall in this class, compared to 21% of Orion sources. This relative 
lack of embedded sources is observed in both Chandra fields. It is also not a property of 
the typically higher mass of our Cyg 0B2 sources compared to those in the deeper low-mass 
COUP sample: filtering both samples on X-ray luminosity (which correlates well with stellar 
mass) and comparing only those sources with Lx > 29.5 erg s~^ (our approximate X-ray 
completeness limit, see Section also shows a deficit of embedded Cyg 0B2 sources. 

While it appears there is a moderate evol ution in stellar X-ray lumin osity with age {Lx 



decreasing with age, but Lx/Lboi increasing, iPreibisch fc Feigelsonll2005l ). and there might 



also be an associated change in spectral hardness as stars age, we expect such effects to be 
very small compared to the gross changes observed in median photon energy as a result of 
the different absorption characteristics of the ONC and Cygnus 0B2 regions. We conclude 
that the relative lack of sources with high median X-ray energies indicates that, unlike the 
ONC, the central region of Cyg 0B2 does not contain an embedded population of recently 
formed stars and that star formation in the vicinity of Cyg 0B2 has declined significantly. 



4.2. The X-ray luminosity function 



Different youn g star forming reg ions, such as Ceph eus 0B3 (iGetman et al.l 120061 ). Ml 7 
jBroos et ahlboOTh . and NGC 6357 JWang et al.ll2007h . have all been observed to have X- 
ray luminosity functions (XLFs) that exhibit similar lognormal shapes. Since the XLF is 
effectively a convolution of the IMF and the X-ray luminosity - mass relation for pre-MS 
stars, it may be used to study the properties of a stellar population without making any 
assumptions about the region or the objects being studied (e.g. by using pre-MS evolution 
models). Furthermore, if we assume that the IMF is invariant in young clusters, at least as 
far as can be measured, the completeness of our sample may be estimated by comparing the 
XLF with that determined from deeper observations of other clusters, such as the COUP 
study which is complete to ~0.1 Mq. 

Figure [5] shows a comparison between the XLF of lightly obscured {{E) < 3.0 keV) 
Cygnus 0B2 sources and those from the well-studied Orion Nebula (we show only the lightly 
obscured populations since Cyg 0B2 has a smaller heavily obscured population compared to 
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Median energy (l^eV) 



Fig. 4. — The distribution of median en ergies of the Cygnus B2 Chandra sources (full line) 
compared to the Orion Nebula Cluster (iGetman et al.ll2005bl . dashed line). 
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Fig. 5. — The extinction-corrected full-band (0.5 - 8.0 keV) X-ray luminosity function for 
all "lightly obscured" {{E) < 3.0 keV) Chandra sources in Cygnus 0B2 (black line) com- 
pared to the equival ent XLF from the Chandra Orion Ultradeep Project (COUP, red line, 
Getman et al.ll2005al) . which has been adjusted to the revised distance of 400 pc for the ONC 
(IMuench et al.ll2008l ). The COUP XLF has been renormalized to fit the high-mass Cyg 0B2 
XLF. All Cyg 0B2 X-ray luminosities are scaled to a distance of 1.45 kpc. XLF Icr error 
bars for the Cyg 0B2 sample are also shown. 
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the ONC, and this will influence the XLF). Error bars for the Cyg 0B2 XLF showing the la 
confidence level were estimated from Monte Carlo simulations of the XLF, where individual 
X-ray luminosities were drawn from a Gaussian distribution with mean equal to the measured 
X-ray luminosity and variances based on the X-ray luminosity error. The COUP XLF has 
been renormalized to fit the Cyg 0B2 XLF at the intermediate- to high-mass end. 

Between log Lx ~ 29.5 - 31.5 erg s~^ the ONC and Cyg 0B2 XLFs are in good 
agreement. At lower luminosities the two XLFs differ, with fewer Cyg 0B2 X-ray sources 
compared to the ONC. This represents our completeness limit, fro m which we determin e a 



mass limit of 1 M0 using the X-ray luminosity to mass correlation in iPreibisch et al. 
Beyond this our completeness drops off rapidly towards lower masses with our sample only 
50% complete at O.75M0 and < 20% complete at O.5M0. Considered separately the two 
fields show similar XLFs down a similar completeness limit, despite an exposure difference 
of a factor of two. This is a product of the Lx - mass relation whi ch is level between 1 



10 Mq and then drops rapidly at lower masses with a high dispersion (IPreibisch et al.ll2005l ). 
Beyond the completeness limit, the XLF falls off quicker in the north-western field than it 
does in the central field due to the different in X-ray exposure and the lack of 



5. Stellar masses and mass functions 

The compilation of a large sample of cluster members with near-IR photometry and 
fitted pre-MS isochrones allows stellar masses to be estimated for a large and relatively 
complete sample. This then allows the mass function of sources in Cyg 0B2 to be determined. 

Individual stellar masses were calculated based on the positions of the sources in the 
(J, J — H) CMD by tracing an extinction vector back to the theoretical MS or pre-MS 
isochrones described in Se ction 13.11 For this purpose we assume a distance to Cyg 0B2 



of 1.45 kpc iHansod ( 120031 ). We will initially assume that all stars in each field were born 
coevally, ignoring any possible age spread, therefore allowing a direct relationship between 
the position of the source in the CMD and its mass and extinction. For this we will use 
the best-fitting ages determined in Section 1371] for DM = 10.8. Because of the form of the 
pre-MS isochrone in the near-IR CMD a number of sources have degenerate solutions for 
their masses. This degeneracy was removed by dereddening these sources using their X-ray 



^Note that a n X-ra y luminosity of log Lx ^29.5 corresponds to a mass of ^0.25 Mq using the conversion 
for Siess et all (l2000l) models presented by IPreibisch et al.l (|2005l ) , but since the X-ray luminosity - mass 
relation ha s a spread of at l east 1 dex, we estimate the actual completeness using the range in luminosity in 
Figure 3 of IPreibisch et al.l |2005l ). 
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spectral-fitted hydroKen colu mn densities and the relation Nh = 2.2 x lO^^Ay cm~^ (e.g. 
Gorensteinlll975l : lRyterlll996l ) and then finding the closest solution on their reddening locus. 
Table [U lists the masses and extinctions for all sources in Cyg 0B2 with near-IR photometry 
determined using this method. 

Masses derived using this method will be dependent on a number of factors. The J- 
band magnitude errors, as shown in Figure 121 are small, and will therefore produce only a 
small error in the derived mass (at J ~ 17, where M = 0.65Mq, the typical J-band error of 
0.024 mags corresponds to an error of only 0.02 Mq). Other photometric "discrepancies", 
such as treating a binary star as a single star or not considering the effects of variability 
will affect the derived stellar mass slightly more. The influence of circumstellar material on 
the stellar masses is reduced by our use of the J — H color as opposed to the H — Kg color 
in deriving masses from near-IR colors. The age of the region, estimated in Section 13. Ij 
will also affect the derived masses, and this is discussed later. Finally, the choice of pre- 
MS isochrones could potentially have a signiflcant effect. Stud ies show that using differen t 
isochrones can produce mass differences of up to a factor two (IHillenbrand fc Whitd l2004l ) . 
though this effect is strongest at low masses. 

For many of the most massive stars in Cyg 0B2 a considerable body of spectroscopic 
observations exists that may be used to derive spectral types and masses. In Paper 1, 
54 sources with spectroscopic observations in the literature were associated with Chandra 
sources of which 51 have full spectral types. We determined zero- ag e masses for a ll these 
sources using the zero-age masses from the published parameter s of iMartins et al.l (120051) 
for 0-type stars and from the collated information presented by Kiminki et al.l (2007) for 
B-type stars (primarily from stellar evolutionary models). 



5.1. The mass function of Cygnus OB2 

The mass functions (MFs) of very young SFRs are often described as initial mass func- 
tions, although this is rarely true. The IMF may be disrupted within as little as ~l-2 Myr 
due to the ejection of members through interactions and the demise of the highest mass 
members. Considering the ages estimated i n Section 13.11 and evidence in the literature for 



ejected massive stars from Cyg 0B2 (e.g. IComeron fc Pasqualil 120071 ). it is important to 



make this distinction and understand its implications. 

Figure [6] shows mass functions for both flelds in Cyg 0B2 and for the combined flelds, 
all of which reach a maximum at ~1 Mq, our estimated completeness limit (Section I4.2p . 
Excluding the bins spanning the range 3-10 Mq that contain A and late B-type stars not 
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Table 1: Derived stellar properties of sources in Cygmis 0B2 with optical or near-IR associ- 
ations. Only the first 10 rows are shown. The complete table is available in the electronic 
edition of the journal. 



No. 


RA 


Dec 


Member Mass 
(Mo) 


Ay 


Source 


1 


20:31:16.03 


41:26:32.2 


Y 


1.9 


10.0 


CMD 


2 


20:31:19.32 


41:29:32.3 


Y 


5.7 


8.9 


CMD 


3 


20:31:20.60 


41:31:51.1 


Y 


5.2 


7.7 


CMD 


4 


20:31:21.15 


41:29:04.1 


Y 


0.66 


6.7 


CMD 


5 


20:31:22.14 


41:28:40.3 


Y 


1.8 


2.2 


CMD 


6 


20:31:22.28 


41:29:29.9 


Y 


1.6 


9.0 


CMD 


7 


20:31:23.32 


41:31:50.7 


Y 


1.4 


3.1 


CMD 


8 


20:31:23.47 


41:27:15.1 


Y 


1.7 


8.1 


CMD 


9 


20:31:23.56 


41:29:48.9 


Y 


1.9 


2.8 


CMD 


10 


20:31:23.63 


41:30:40.7 


Y 


1.3 


5.4 


CMD 



Notes. Col. (1): X-ray source number. Cols. (2)-(3): Position of optical or near-IR source 
that X-ray source was matched to. Col. (4): Association membership information (Y/N). 
Col. (5): Stellar mass. Col. (6): Extinction toward source. Col. (7): Method by which mass 
and extinction were determined. For sources with spectral information in the literature a 
reference is provided, while sources where masses were determined from the CMD are listed 
as "CMD". 
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expected to emit X-rays (see Section [2] and our discussion of the completeness of our sample) 
the mass functions show a clear power law slope, steepening at higher masses. Excluding the 
most massive stars we find a power-law slope of F = —1.0 9 ± 0.13 i n the combined fields, in 
agreement with the value of F = — 1.3 ± 0.3 estimated by iKroupal (120011 ) from accumulated 
observations, and only slightly shallower than the canonical Salpeter value. When the two 
fields are considered separately we find negligibly different slopes of F = —1.08 ± 0.15 in the 
central field and F = —1.09 ± 0.10 in the north-western field. When considering only the 
0-type stars we find a much steeper slope with F = —2.72 ± 0.52. 

Previous measurements o f the MF§| in Cyg 0B2 have mostly been derived from spec- 
troscopy of the massive stars. iMassey &: Thompson! (jl99ll ) measured F = —1.0 ± 0.3 for all 
stars more massive than 9 Mq and F = —0.9 ± 0.5 using only the highest quality photom- 
etry and stars more massive than 15 Mq. Compared to our high-mass MP this is a much 
shallower slop e, although the masse s used were derived from older models for massive stars. 
More recently, iKiminki et al.l (120071 ) derived a slope of F = —2.2 ±0.1 for spectroscopically 
studied stars more massive than 15 M0 across a wide area of Cyg 0B2, deriving masses from 
updated stella r evolu tionary models and the revised O star stellar parameters presented by 
Martins et al.l (120051 ). Our high- mass MP is in good agreement with their value, mainly 
because of our use of similar models and stellar parameters. However, we have shown here 
that the MP at lower masses is not as steep and is closer to the canonical value. We suggest 
that the steep high-mass slope of the MP is not intrinsic but is a product of the older stellar 
populations in the region having lost their most massive members, thereby steepening the 
mass function. 



The only previous measurement of the mass function at lower masses is due to lKnodlseder 

(12000! ) who found F = —1.6 ±0.1 down to 2 Mq from a photometric analysis of 2MASS obser- 



vations in Cyg 0B2. However, his use of main-sequence magnitudes and colors for a pre-MS 
stellar population puts the derived masses into question. Our measurement is therefore the 
firs t deep na e asure ment of the pre-MS mass function in Cygnus 0B2 and its agreement with 
the iKroupal (I2OOII ) value, supporting the view that the IMP is invariant. 



•^The majority of previous mass function measurements for Cyg 0B2 in the literature refer to the IMF, 
most hkely because of the previously assumed age of ~2 Myrs. A full understanding of the age and star 
formation history of the region is therefore clearly important for correctly interpreting the measurements of 
the mass function. 
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5.2. Effect of age uncertainties on mass function slope 

The uncertainty in our estimation of tlie age of Cyg 0B2 fias repercussions for stellar 
masses derived from fitting pre-MS isochrones to the near-IR CMD. We also cannot assume 
that all the stars within a single field were formed coevally (see the discussion in Section |6TT]) . 
To test the impact on the MF slope of this age uncertainty we altered the estimated age and 
regenerated the individual stellar masses and the MF. We find that when the estimated age 
is varied by ±1 Myr the derived MF slope changes by less than 10%, with the MF in the 
north-western field varying much less than in the central field. Even varying the age over a 
large range (by ±2 Myr) only causes a MF slope variation of ~20% (for the central field) 
and < 10% in the north-western field. Age variations larger than this cause much greater 
variations in the mass function slope, particularly in the central field, but are well outside 
the 95% confidence intervals. 

We also studied how the MF would vary if all the stars were not formed coevally but 
formed over a range of ages (an age spread that might explain the wide spread in the CMD, 
see discussion in Section 16.11) . Since the spread in the CMD could be caused by both a 
spread in extinction and a spread in ages we did this by selecting ages for each sources 
that would minimize the resulting spread in extinction. This is the opposite process to that 
used above where a single age was used, thereby minimizing the spread in ages (the true 
situation is likely to be somewhere between these two extremes). The resulting MF has a 
slope of r = —1.03 ± 0.11 in the combined fields, F = —0.98 ± 0.13 in the central field, and 
r = —1.05 ± 0.06 in the north-western field. These small variations in the slope of the MF 
strengthen its validity despite the uncertainty over the age of the region. 

6. Discussion 

In this section we discuss the results obtained here in the context of our current un- 
derstanding of Cygnus 0B2, its size, age, and star formation history. We also compare the 
derived properties of the region with those of the small but growing number of massive SFRs 
that have been studied. 



6.1. The age of Cygnus OB2 



The main constraints we hav e on th e age of Cyg 0B2 are derived from observations 
of its most massive stars. iHansonl (120031 ) estimated an age of 2 ± 1 Myrs to accommodate 
the observed mix of main sequence and evolved OB stars: very young sources such as #A37 
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from IComeron et al.l (|2002[ ). which iHansonl (120031 ) classified as an 05V((f)), cannot be older 



than 1 Myr (based on evolutiona ry tracks); while supergiants such as Cyg 0B2 #9, type 
03If (IMassey fc Thompson! 1 1 99 ll ). are believed to be at least 3 Myrs old. The presence of 



a main sequence extending from type 06, which also cannot be older than 3 Mys, supports 
this age estimate. 

Th ere has been growing evidence for a large numbe r of evolved supe rgiant stars in the 
region (IMassey fc Thompson! !l99l! : !Comer6n et al.! !2002! ) . !Hanson! (12003! ) noted that as the 
search for OB stars in the region was extended to larger radii the fraction of evolved massive 
stars increased, leading her to suggest that our view of Cyg 0B2 is b ecoming contarninated 



by non-members. Many of the new evolved massive stars identified b v!Comer6n et al 



were found to the south of the original association in a region where !Drew et al.l (!2008!) have 



20021 ) 



Drew et al. 



since uncovered a population of A-type dwarf stars with an age of 5 - 7 Myr 
(12008! ) suggest that this region may represent a larger and older parent cluster. 

In Section !3?T! we derived ages of 3.511i'q^ and 5.25^} o Myrs for the central and north- 
western fields, both with large spreads in the CMD that cannot solely be explained by 
photometric errors. Such wide spreads in CM Ds have previously be en interpreted as due to 



wide age spreads in association members (e.g. !Palla &: Stahlerl!2000[ ). However other authors 



have suggested that the observed CMD spreads are due to a mixture of photo metric errors, 



variability in pr e-MS stars (up to several mag nitudes for classical T-Taur i stars. !Herbst et al. 



1994! ). binarity ( IPreibisch fc Zinnecker!!l999! ). episodic accretion phases (IBaraffe et al. 



20091 



or a spread in line-of-sight distances or extinctions. Many of these factors a re not well under- 
stood, and in the case of Cyg 0B2 particularly the spread in extinction. !Sale et al.! (120091 ) 
measured the rise of extinction with distance towards Cyg 0B2 and found a substantial 
increase in extinction at ~1.5 kpc, which they estimate could be confined to a depth of 
< 100 pc. If a significant spread of extinction did exist across Cyg 0B2, it could, combined 
with the effects of binarity and pre-MS variability, explain a large proportion of the spread 
in the CMD. Furthermore, given the uncertainties on the two ages derived in this work, the 
age difference between the two fields is not significant in itself and could fit with a picture 
of a single burst of star formation in the region. 

However, combining the ages derived here with the existing evidence for very young 
(~1 Myr) and slightly older (~5 - 7 Myr) populations in Cyg 0B2 suggests a picture of a 
region with a large spread in ages. In such a picture the central field, for which we derived 
an age of 3.5t[f Myr, merges with the young OB-star dominated association, while the 
north-western field is biased to a greater age. We cannot yet know to what extent there 
are two distinct generations of stars present, or whether there is a real age spread (with the 
youngest stars co- located with the OB stars). 
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6.2. The size and mass of Cygnus OB2 



As shown in Section [2l the catalogue presented in Paper 1 is complete to 1 Mq within 
the caveats of X-ray emitting stars. We can therefore combine the MF determined here 
and the stellar density within our completeness limits to estimate the total size and mass 
of the entire Cyg 0B2 association. For the purposes of this exercise we will consider the 
Cygnus 0B2 association to consist of all sources young enough to be detected in X-rays (i.e. 
< 100 Myrs). While this may include multiple generations of star formation that may have 
occurred on different spatial scales, until the different populations can be separated this is 
the only option available to us. Also, because our observations only sample two small fields 
in the center of the association, we must combine our results with previous estimates of the 
total spatial extent of the SFR available in the literature. 



Knodlsederl (120001 ) analysed the morphology of Cyg 0B2 by fitting a King profile ((King 
19621 ) to 2MASS star count data, estimating it to be a spherically symmetric association 
of ~ 2° diameter. This analysis necessarily relied heavily on performing a background 
subtraction using the stellar content of a nearby field, despite the inherent difficulties involved 
in doing this along complicated sightlines in the Galactic plane. An improved morphological 
analysis may be made using the radial stellar density of Chandra sources in the direction 
of the north-western field where a radial length of 30' is possible using the two fields (see 
Figure 1 in Paper 1). 



We find that the stellar density drops off much more quickly than estimated by lKnodlseder 



(I2OOOI ). with approximately one third as many stars detected by us in the 1 - 2.5 Mq mass 
range (where our d ata remain c omple te) at distances of >15', compared with the King pro- 
file determined by iKnodlsederl (j2000( ) . Using our radial density profile and assuming that 
Cyg 0B2 is spherically symmetric we have estimated a new King profile for Cyg 0B2. Using 
our peak stellar density and the mass function slope e stimated above, extrapolated down 
to 0.01 Mq using the multi-stage power-law IMF from iKroupal ( 1200 ll ) we calculate a total 
stellar content of 5.7 x 10*^ stars and a total mass of 2.8 x 10^ Mq. However it is clear from 



the spatial distribution of A-type dwarfs in Cyg 0B2 presented by iDrew et al.l (120081 ) that 
the structure of Cyg 0B2 is far from spherical. The authors note a distinct lack of A-type 
dwarfs to the east of Cyg 0B2 and a pronounced concentration towards the southern part 
of the cluster that may increase the size of the association. We therefore estimate a total 
stellar mass of (3 ± 1) x lO'* Mq for the entire association. 

Using this estimated density profile and our derived MF we can estimate the expected 
number of stars of different masses in Cyg 0B2. A total mass of 3 x 10^ Mf;^ implies a total 
of ~1200 OB stars, approximately half that predicted by IKnodlsederl ( I2OOOI ). Our current 
MF predicts a total number of 0-type stars of ~75, only slightly larger than the currently 



- 23 - 



known number of 65 0-type stars in Cyg 0B2 (jMassev fc Thorapsoru Il99ll : iHansonl 12003 



Kiminki et al.l 120071 : IComeron et al.l l2008l : iNegueruela et al.l |2008| ) . This may suggest that 



there are more 0-type stars to be discovered in Cyg 0B2. If we assume that the high-mass 
MF slope has been influenc ed by the dem ise of the most massive stars in the region and that 
it was originally F = —1.3 (jMasseylll998l ). we estimate the number of stars that have evolved 
to their end states to be ~12. There are no known sup e rnova rem nants in Cyg 0B2, although 
there is some evidence for stellar remnants: iBednarekl (120031) and 
two pulsars within Cyg 0B2, with ages 



2 X 10"^ yr (iCamilo et a 



Abdo et all (l2009f ) identify 



20091 ) , that likely represent 



recently expired massive stars, while Comeron &: Pasqualil ( 2007 ) identified an 04If runaway 
star that appears to have been ejected from Cyg 0B2 ~1.6 Myr ago, possibly following the 
supernova explosion of a more massive binary companion. 



6.3. Comparison with other massive star forming regions 



Table [2] lists some of the properties of Cyg 0B2 estimated in this work along with those 
of a number of other massive star forming regions. With the exception of the ONC the 
majority of these regions have not been studied to as low a stellar mass, with properties 
determined in different ways by different authors (the number of 0-type stars is particularly 
poorly known for many of these regions). Despite this some comparisons may be made. 
Cyg 0B2 is certainly a very high-mass cluster, comparable to the Arches and Quintuplet 
clusters near the Galactic center or the giant H ll regions Westerlund 2 and NGC 3603 in 
the Carina spiral arm. Only regions such as Westerlund 1 cluster, W49A, and R136 in the 
Large Magellanic Cloud appear to be more massive. 

Despite their similar masses, these SFRs show significant differences in their stellar 
densities, ranging from ~ 10^ Mq pc~^ for Cyg 0B2 and the Quintuplet cluster to ~ 10^ Mq 
pc~^ for NGC 3603 and the Arches cluster. This may represent an inherent difference between 
two varieties of massive star forming region: that of the compact, massive cluster with a very 
high density core on one hand and the lower density and spatially more extended region on 
the other. This could be an evolutionary effect since Cyg 0B2 and the Quintuplet cluster 
are both older than NGC 3603 and the Arches cluster, alt hough the age differenc es are much 
shorter than the time periods necessary for such a change (IGieles fc Bastianll2008l ). It is more 
likely that these differences are real and are related to how the molecular cloud collapsed: 
low density cores suggest formation by triggering or the interaction o f gas motions, while 
gravitational collapse typically results in higher densities (ISmithll2008l ). 



Homeier fc Alved (120051 ) note that W49A, which also has a low core density, has a very 



sub clustered structure, a picture that could apply to the wider Cygnus X giant molecu- 
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Table 2: Comparison of the properties of Cyg 0B2 with other regions. The number of 0-type 
stars for Cyg 0B2 are only those spectroscopically confirmed, while for other regions they 
are generally extrapolated from photometry. 



Star forming 


Age 


Mtotal 


Core density 


Core radius 


stars 


region 


(Myrs) 


{Me) 


(Mo pc-3) 


(pc) 




ONC " 


0.3-1.0 


1.9 X 10^ 


2 X 10^ 


0.16-0.21 


1 


Westerlund 2 


2.0 ±0.3 


7.0 X 10^ 




0.002-0.2 


> 24 


Quintuplet Cluster'^ 


4.0 ±2.0 


6-10 xlO^ 


2.5 X 10^ 


1.0 




NGC 3603 ^ 


1.0 ±0.35 


1.0 X 10^ 


1 X 10^ 
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lar cloud (GMC) wit h subclusters such as Cyg 0B2 Jsica et allboosl ) and Cyg 0B9 (e.g. 
Uyaniker et aLll200ll ). This structure supports the hierarchical picture of star formation 
where young clusters and associations represen t the dense inner region s of a hierarchy emerg- 
ing from the GMC in which they formed (e.g. lElmegreen et al.ll2006l ). 



Conclusions 



We have analysed the integrated stellar properties of young stars in the massive SFR 
Cygnus 0B2 using the recently published catalogue of Chandra X-ray point sources detected 
in the region. We find that our sample is complete down to ~1 Mq, excluding A and B-type 
stars that are not expected to emit X-rays, making this one of the deepest surveys of the 
region. Ages for sources in the two fields studied were estimated from pre-MS isochrone 
fits to the near- IR CMD and found to be 3.5l?;o^ and ^.2ht\l Myrs, both with significant 
spreads around the isochrones that cannot solely be explained by phot ometric errors. Given 
that t he most mass ive stars in the region have an age of 2 ± 1 Myrs (IMassey fc Thompson 
199ll : iHansonl |2003| ) . we su ggest that this supports recent eyidence for an older generation 



of stars in the region (e.g. iDrew et al.l 120081 : IComeron et al.l |2008[ ) that either represent a 



significant age spread or multiple epochs of star formation. Additionally, the lack of sizable 
population of highly embedded sources suggests that star formation in the region studied 
has declined significantly. 

Assuming that a significant fraction of sources in Cyg 0B2 are older than the pre- 
viously accepted age of 2 Myrs, many of the unique properties of the region may be ex- 
plained. As noted by lDrew et al.l ( l2008l ). the distinct lack of a bright Hll region surrounding 
Cyg 0B2 is perfectly reasonable if the clearing timescale has been greater than 5 Myrs. 
Furthermore, we propos e that the low fraction of sources with i^'-band excesses found by 
Albacete Colombo et al.l ( 120071 . and confirmed from our larger sample) is less abnormal if a 
significant fraction of sources are older than previously believed. 

We measure the stellar mass functions for stars within our completeness limits and 
find a slope of F = —1.09 i fc 0.13, in good agreement with the potentially universal value 



estimated by iKroupal ( 12002| ) . This is the deepest and most accurate measurement of the 
mass fu nction in Cyg 0B 2, building o n previously measurements that have found either 
steeper ( iKnodlsederl 120001 ) or shallower ( iMassey fc Thompsonlll99ll) slopes . A st eepening of 
the mass function is observed at high masses, as found by iKiminki et al.l (120071 ). which we 
attribute to the demise of the highest mass members of the older generations in the region. 



An estimation of the total size of the region suggests that the stellar density drops 
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off faster tlian estimated by iKnodlsederl (120001 ) . Tliis, combined witli our mass function 



estimate, reveals a total mass of ~ 3 x 10^ Mq. This makes Cyg 0B2 similar in mass 
to many of the recently discovered massive clusters in our Galaxy such as the Arches and 
Quintuplet clusters. However, we note that comparing other properties of these massive star 
forming regions, such as their central densities, reveals differences that may hint at different 
formation mechanisms. 

Past studies of Cyg 0B2 have clearly been complicated by an uncertain star formation 
history and a large spatial extent in a crowded region of the Galactic Plane. These are likely 
to be inherent properties of massive star forming regions in galactic spiral arms where star 
formation continues for longer periods and across larger areas, propagated through entire 
GMCs by continued rejuvenation and feedback triggering. These findings support the need 
to fully understand the full star formation history of a region before if its observed properties 
are to be correctly interpreted. Cygnus 0B2 is clearly a very important region to study, not 
only because of its proximity and size, but also because of the insights it will give us into 
galactic-scale star formation and the evolution of GMCs. 
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M / Mg M / Mg 

Fig. 6. — Mass functions (MF) for all sources in Cygnus 0B2 with stellar masses derived 
from near-IR photometry and existing spectroscopic observations. Left: MF for both fields 
combined; right: MF for the central (full line) and north-western (dashed line) fields, 
fits to the MFs within the completeness limits (M > IMq and excluding early A and late 
B-type stars) are shown as red lines, with slopes of F = —1.09 ± 0.13 (combined field), 
F = -1.08 ± 0.15 (central field), and F = -1.09 ± 0.10 (north-western field). 



